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Abstract 

Recently, osmium diboride (OSB2) has attracted considerable attention as an incompressible and 
hard material. We investigate the structural property, elastic constant, and electronic structure of 
orthorhombic OsB2 by the first-principles total energy calculations. The calculations are performed 
within the density functional framework using the projector augmented plane wave method. The 
structural property and bulk modulus of OsB2 compare well with experimental data. The nine 
independent elastic constants of orthorhombic OsB2 at zero-pressure have also been calculated by 
symmetry-general least-squares extraction method. We have analyzed the mechanical stability of 
orthorhombic OsB2 in term of the calculated elastic constants. A detailed study of the electronic 
structure and the charge-density redistribution reveals the features of strong covalent B-B and Os- 
B bondings in orthorhombic OsB2- The orbital hybridization and the characteristics of bonding 
orbitals in OsB2 are identified. Orthorhombic OsB2 exhibits a metallic character and the states at 
Fermi level mainly come from the d orbital of Os atoms. 
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I. INTRODUCTION 



To design and explore ultrahard materials is of great scientific and practical interest 
in both experimental and theoretical studies]]], §, £|, ^ 0, because the class of hard 
materials has a wide variety of important technological applications ranged from cutting and 
polishing tools to wear-resistant coating ||, [/] . Stiffness reflects the elastic compressibility 
of materials, which is related to the bulk modulus, shear modulus, and the Young's modulus. 
Hardness describes the resistance of materials to plastic deformation 0. It is recognized 
that there is a good correlation between shear modulus and hardness [Q, [5], |8| and a three- 
dimensional network composed of short , strong bonds is critical for hardness of material || . 
Many synthetic and theoretical efforts are made to design new ultrahard materials in past 
decades and two main approaches have been summarized || as the following: a) combination 
of light elements (such as boron, nitrogen, and/or oxygen) to form short covalent bonds and 
b) compounding of elements with very high densities of valence electrons to ensure that the 
materials resist being squeezed together. Cubic BC2N, an example of the first approach, was 
synthesized as a dense ternary phase in the B— C— N system [7] and found to be an extreme 
hardness second only to diamond [0, [|. For the second approach, one usually uses transition 
metals having a high bulk modulus but low hardness to combine with small, covalent bond- 



forming atoms 0. For example, highly incompressible phase of RUO2 ]T0[ owes to a strong 
covalent bonding between ruthenium d states and oxygen p states [[□]] . 

Other representative example for the second approach is the compounds related osmium 
(Os) element. The pure Os metal was measured to have a high bulk modulus of 462 GPa 



||12|| , 411±6 GPa [fL3|j , and 395±15GPa [|1J] by different experimental methods, but found to 
have a low hardness of only 400 kg/mm 2 ||15|1 . These features of could be understand that || 
the high bulk modulus of Os mainly owes to the high valence electron density (0.572 e/A 3 ) 
and the low hardness of Os metal with hexagonal close-packed crystal structure is related to 
the metallic boning character. To improve the hardness of osmium, incorporation of small 
and covalent main group elements (such as boron, carbon, nitrogen, and/or oxygen) into 
osmium metal has been proposed in experimental and theoretical studies || [U], [If], [17) . 
Os02 is expected to have a covalent bonding and a high bulk modulus fllf , analogous to 
Ru0 2 . Although OsC and OsN compounds are not yet synthesized successfully, recent ab 
initio calculations predict that |16|] OsC in hexagonal WC structure has a bulk modulus 



2 



of 396 GPa and is a promising superhard material. For one of binary compounds in Os-B 



system, OsB2, recent experiment measures show that [|Tj]] OsB2 is an ultra-incompressible 
and hard material with bulk modulus of 365-395 GPa and hardness of > 2000 kg/mm 2 . 

In order to investigate the mechanical properties of OsB 2 and understand its high stiff- 
ness and hardness, we have evaluated the nine independent elastic constants including bulk 
modulus and shear modulus of OsB 2 using ab initio pseudopotential calculations. Our re- 
sults show that the behavior of this compound can be understood on a fundamental level in 
terms of their electronic band structure. The unusual hardness originates from the strong 
covalent bonding character between the p bonding states of boron atoms and the d orbitals 
of osimum atom, and the strong covalent bonding between boron and boron atoms. The 
ultra-incompressibility exhibits anisotropy (C 33 > C\\ > C22), which is due to the different 
bonding character of interlayer along c, a, and 6-directions of crystal. 

II. METHOD 

Because crystals deform almost in a linear elastic manner at small strains, the increment 
of strain energy density of a homogeneous and elastically deformed crystal is given by 

dE = Uijdeij (1) 

where are the elements of stress tensor and are the elements of strain tensor. Indices 
i, j, k and I run from 1 to 3. Elastic constants are defined as 

_ dan 1 d 2 E 

Cijkl = -7; — 77n K OijtTkl — Lsijkl — 0ijO~kl \£) 

where V and 5 and are the volume of the unstrained crystal and Kronecker delta function, 
respectively. At zero hydrostatic pressure c^-jy = C^ki- The internal energy E(V, {ejj})of a 
crystal under a general strain can be expressed by expanding the internal energy E(V) 
of the deformed crystal with respect to the strain tensor, as 

3 3 

E(V,{e ij }) = E(V,0) + Vj2a ij e ij + - ^ C ljkl e l3 e kl + . . . (3) 

i,j=l i,j,k,l=l 

where E(V, 0) is the corresponding ground state energy. If one uses Voigt notation (xx — > 1 
, yy — ► 2, zz — > 3, yz — > 4, xz — > 5 and xy — > 6) for the strain tensors subscripts (ij, kl) in 
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the the above expression, Eq. 



can be rewritten as 

6 
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E(V, {e,}) = E(V, 0) + Vj2 + y J2 C ^ + ' ' • 



(4) 
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with the strain tensor 
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Therefore, one can apply small elastic strains and calculate the change of energy or stress 
to obtain elastic constants. Direct calculation of elastic constants are possible from first- 
principles total energy calculations and symmetry-general least-squares extraction method 
[Tj|. The independent elastic constants for crystals can be determined by selectively im- 
posing strains either individually or in combination along specific crystallographic direc- 

|. Table | illustrates the deformations employed in determining 



tions PL 20, 21, 22 



independent elastic constants of orthorhombic crystal fl9], p0| . 

First-principles total energy calculations on OsB 2 were performed using a plane-wave 
pseudopotential method as implemented in the VASP code [p4 , p5fl . The exchange- 
correlation energy functional was treated with the local density approximation (LDA) given 
by Ceperley and Alder p6| as parameterized by Perdew and Zunger f2~7 |. Electron-ion inter- 
action was represented by the projector augmented wave (PAW) method |28], and wave 
functions were expanded by the plane waves up to an energy cutoff of 500 eV. Brillouin- 
zone integrations were approximated using the special fc-point sampling of Monhkorst-Pack 
scheme |50| with a 9x13x9 grid. The total energy convergence test showed that convergency 
to within 1 meV/atom was achieved with the above calculation parameters. The bulk unit 
cell lattice vectors and atomic coordinates of orthorhombic OsB 2 were relaxed at a series of 
fixed volumes by optimizing both the forces and stresses. The residual force and stress in 
the equilibrium geometry are of the order of 0.01 eV/A and 10~ 3 GPa, respectively. The 
obtained energies were fitted with the third-order Birch- Murnaghan equation of state (EOS) 



1 3jJ to give the equilibrium volume and the minimum energy. The final calculated cell pa- 
rameters are given in Table [H| which are in good agreement with the available experimental 
values. 
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III. RESULTS AND DISCUSSIONS 



A. Structural properties 

The single-crystal of OsB 2 was found to have the orthorhombic symmetry with the space 
group Pmmn (No. 59) and lattice constants a = 4.6832 A, b — 2.8717 A and c = 4.0761 
A, using the ploycrystalline x-ray diffraction techniques. ||32|| In the orthorhombic unit cell 
of OsB2, the Os atoms were determined at the Wyckoff position of 2a (j, |, z) and 4/ 
j, v) for B atoms P2] . The crystal structure of OsB 2 is illustrated in Fig.[TJ, especially, 
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the layer structure repeated along c-(or z-)direction with Os-B-B-Os sandwich structure 
is also shown. The orthorhombic unit cell of OsB2 contains two formula unit (f.u.). The 
calculated lattice constant are a = 4.6444 A, b = 2.8505 A c = 4.0464 A, which are in 



good agreement with the available experiment values |32|] and previous spin polarized LDA 
(namely, LSDA) calculation results (a= 4.6433 A, b = 2.8467 A c = 4.0432 A) @. We 
should pointed out that we have also performed the spin polarized LDA calculation on OsB2 
and found that there is no difference between the results of LDA and LSDA. The calculated 
volume (53.57 A 3 /f.u.) is 2.2% smaller than experimental value, which is typical for the 
LDA approximation to DFT. The calculated atomic coordinate parameters are z = 0.1538 
for Os atoms, u = 0.0588 and v = 0.6378 for B atoms. We have also calculated the equation 
of state (the pressure versus volume) by taking the volume derivative of the fitted total 
energy obtained from our LDA calculation, the result along the available experimental data 
is shown in Fig. 0. The bulk modulus and its pressure derivative for OsB2, obtained by 
fitting the total energy curve to third-order Birch- Murnaghan equation of state (EOS) fl3~Tl . 
are 336.1 GPa and 4.27 in our LDA calculations. This bulk modulus is 7.9% ~ 15% smaller 
than the experimental value (365-395 GPa) |l~7| . Note that the corresponding derivative of 



bulk modulus to pressure in experimental studies [17]] is fixed in the range between 4 and 1.4. 
In order to further asses the effects of the exchange-correlation approximation, we have also 
studied structural properties of OsB2 using the generalized gradient approximation (GGA) 



34J. GGA consistently yields a larger volume than experiment, while LDA consistently 
gives a smaller volume than GGA (see Table 0). Since the GGA calculation gives a smaller 
bulk modulus (303.45 GPa) of OsB 2 than the experimental data, in the following sections 
we will mainly present and discuss the results obtained by LDA calculations. 
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B. Elastic constants 



The nine independent elastic constants for 0sB2 were determined by imposing three dif- 
ferent strains given in Table [I] on the equilibrium lattice of orthorhombic unit cell and fitting 
the dependence of the resulting change in energy on the strain. Up to seven deformations, 7 
from -0.009 to 0.009 in steps of 0.003, were applied to determine each elastic constant. The 
ground state energies of OsB 2 unit cell were determined for all these deformation strains. 
The E(j) curves are well fitted by the third-order polynomials in 7 (Eq.([|)), as can be 
seen from the small standard errors in the calculated Cy. The elastic constants were then 
calculated using the energy relations listed in Table |. During the calculations on all elastic 
constants of orthorhombic OsB 2 , relaxation of the coordinate of ions in the strained lat- 
tice was carried out for achieving accurate results. The nine independent elastic constants 
for orthorhombic OsB 2 are listed in Table |HI| , which experimental data on elasticity is not 
available, except for bulk modulus. The nine independent elastic constants of orthorhombic 
OsB 2 are C n = 597.0 GPa, C l2 = 198.1 GPa, C 13 = 206.1 GPa, C 22 = 581.2 GPa, C 23 = 
142.6 GPa, C 33 = 825.0 GPa, C 44 = 70.1 GPa, C 55 = 212.0 GPa, and C m = 201.3 GPa in 
our calculations. 



Under the Voigt approximation |19|, |35[, the elastic constants for polycrystalline aggre- 
gates are described by Voigt shear modulus (Gy) and Voigt bulk modulus (By). For the 
orthorthombic lattice, Gy and By are expressed in the following equations ||19||: 



Bv = ~(Cn + C 22 + C 33 + 2C 12 + 2C 13 + 2C 23 ). (7) 
y 

The calculated Voigt shear constant of OsB 2 is about 193.8 GPa. The Voigt bulk modulus 
calculated from the theoretical values of the elastic constants is 344.1 GPa. It agrees well 
with the one extracted from the fit to the third-order Birch-Murnaghan equation of state, 
336.1 GPa. Since the Young's modulus (E) and Poisson's ratio (v) for an isotropic materials 



can be obtained from the following relations [19 



e = WTg (8 » 

3B-2G 

'~ 2(35 + ^ 



The calculated Young's modulus and Poisson's ratio of OsB2 are 489.5 GPa and 0.263, 
respectively. 

The requirement of mechanical stability in a orthorhombic crystal leads to the following 
restrictions on the elastic constants: PO, 



(C xl + C 22 - 2C 12 ) > 0, 
(Cu + C33 - 2C 13 ) > 0, 
(C 22 + C 33 - 2C 23 ) > 0, 
Cii>0,C 22 >0,C 3 3>0, (10) 

C44 > 0, C55 > 0, C*66 > 0, 

(C u + C 22 + C 33 + 2C 12 + 2C 13 + 2C 23 ) > 0. 

Since B is a weighted average of Cu and CV,- (indices i and j running from 1 to 3), these 
conditions also lead to a restriction on the magnitude of B . Namely, B Q is required to be 
intermediate in value between \{C\ 2 + Ci 3 + C 23 ) and \{Cxi + C 22 + C 33 ) 0, |36| , 



\{C l2 + C 13 + C 23 ) <B < i(C n + C 22 + C 33 ) (11) 

The elastic constants of orthorhombic OsB 2 in Table [111] obey these stability conditions 
(Eq. (|I0|) and (|TT]) ) . In particular, C\ 2 is smaller than the average of C\\ and C 22 , Ci 3 is 
smaller than the average of of C\\ and C33, C 2 3 is smaller than the average of C 22 and C33, 
and the bulk modulus is smaller than the average of Cu, C 22 , and C33 but larger than the 
average of Ci 2 , C13, and C 23 . 



C. Electronic structures and analysis of chemical bonding 

The energy band structure and density of states (DOS) (including total, site and angular 
momentum decomposed DOS) of orthorhombic OsB 2 are shown in Fig. |2| and ^ respectively. 
It evidently exhibits a metallic character. The valence bands in Fig. ||] have a width of about 
15 eV, and are split into two disjointed groups, respectively. The lower group has a width 
of about 5 eV and the upper one up to the Fermi level (Ep) has a width of about 10 eV. 
Combining from the DOS figures (|3|), the lower group in the energy range between -15 eV 
and -10 eV below the Fermi level {Ep) is largely contributed by B-2s levels and small amount 
of Os-6s, the upper one in the energy range from -12 eV to Ep is mainly composed by Os-5<i 
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orbitals and B-2p. It suggests that the strong hybridization between B-2p and Os-5d orbitals 
occur below the Fermi level. Additionally, d orbitals of Os atom are contributed most to 
the total DOS at E F . These findings are consistent with the previous studies on electronic 
structure of OsB 2 [KBI, 



The large values and anisotropy of elastic constants of OsB 2 can be understood by ex- 
amining the nature of chemical bonding between the atoms. In order to understand the 
bonding between the Os-B, Os-Os, and B-B atoms, the lines charge density along nearest 
neighbor Os— B, Os— Os, B— B atoms are illustrated in Fig. [|. It is clearly seen that charge 
highly accumulates around Os atom, and Os-B, Os-Os and B-B bondings exhibit covalent, 
metallic, covalent characters, respectively. It also supported by above discussion that the d 
states of Os atoms and p states of B atoms have a strong hybridization seen in density of 
states. In order to further reveal the topology of Os-Os, Os-B and the B-B bondings, we 
illustrate the contour plots of electron density in different crystallographic planes. Fig. |6| 
illustrates the two-dimensional electron charge density distribution in the planes with nor- 
mal vector along <100>, <010> and <001>, respectively. From Fig. ^(a) and (d), in 
which planes with the normal vector along <100> and central point at Osl atom (|, |, 
z), and with the normal vector along <001> and central point at Osl atom, respectively, 
it can be clearly seen that there is no strong bonding between Os— Os atoms as there is 
no accumulation of charges in excess of the background charge density of the plane. It 
appears that the Os— Os bonding is predominantly metallic in OsB 2 . Fig. ^|(e) and (f), in 
which planes with the normal vector along <001> and central point at Bl atom (u, |, v), 
and through Bl, B2, B3 and B4 atoms, respectively, It also can be clearly seen that there 
is strong covalent bonding between B-B atoms. The nature of Os-B and B-B bonding is 
illustrated in Fig. P(c) of the two-dimensional contour electron density map with normal 
vector along <010> and central point at Osl atom (7, |, z). It can be seen that there is 
a charge accumulation of electron charge between the Os-B atoms, which is consistent with 
the idea of B-2p and Os-5d hybridization noted in Fig. [|. This also clearly confirms both 
B-B and Os-B bondings are covalent character. Because the orthorhombic OsB 2 is stacked 
by Os-B-B-Os atomic layer structures in the c (or z) direction of crystal, these mean that 
the atom bondings of interlayer are strong covalent characters of Os-B and B-B. This also 
explains the high value of elastic constant C33 (825.0 GPa) and a relatively lower value of 
elastic constant C n (597.0 GPa) and C 22 (581.2 GPa). 
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IV. CONCLUSIONS 



The structural parameters, elastic constants, and electronic structures of OsB2 have been 
calculated and analyzed. Our calculated structural parameters including the lattice con- 
stants and internal parameters are in good agreement with the experiment data. The LDA 
calculations yielded the bulk modulus (336.1 GPa) of OsB2 , which is more close to the ex- 
perimental data than that of GGA calculations (303.45 GPa). The LDA calculations also 
predicted that the independent elastic constants of OsB 2 are Cxi = 597.0 GPa, C\% = 198.1 
GPa, C 13 = 206.1 GPa, C 22 = 581.2 GPa, C 23 = 142.6 GPa, C 33 = 825.0 GPa, C 44 = 70.1 
GPa, C 55 = 212.0 GPa, and C 66 = 201.3 GPa, respectively. The incompressibility of OsB 2 
exhibits anisotropy (C 33 > Cn > C 22 ) and the c-direction of crystal is the least compress- 
ible due to the difference of atom bonding character of interlayer along different crystalline 
planes. The detailed study of the electronic structure and charge density redistribution 
reveals the strong covalent Os-B atoms bonding and B-B atoms bonding play a important 
role in the incompressibility and hardness of OsB 2 . 
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TABLE I: Parametrization of the strains used to calculate the elastic constants of orthorhombic 
OsB 2 . The energy expressions reported on the third column refer to Eq. ||. 



Strain 


Parameters (unlisted e, : 


= 0) 




AE/V in C>(7 2 ) 


e 1 


ei = 7 






he*? 


e 2 


e 2 = 7 








e 3 


e 3 = 7 








e 4 


ei = 27, e 2 = -7, e 3 = 


-7 


i(4Cu - 4C 12 


- 4C 13 + C 22 + 2C 23 + C 33 ) 7 2 


e 5 


ei = -7, e 2 = 27, e 3 = 


-7 


UCu - 4C 12 - 


1- 2Ci 3 + 4C 22 - 4C 23 + C 33 )7 2 


e 6 


ei = ~7) e 2 = ~7; e 3 = 


2 7 


^(Cn + 2Ci 2 


- 4Ci 3 + C 22 - 4C 23 + 4C 33 )7 2 


e 7 


e 4 = 7 






|C 44 7 2 


e 8 


es = 7 








e 9 


e 6 = 7 






lC 667 2 



TABLE II: Calculated lattice parameters (a, b and c, in A), volumes of unit cell (V, in A 3 ), and 
bulk modulus (B, in GPa) of orthorhombic OsB 2 , along a comparison with other theoretical work 
and available experimental data. 



Property 


This work 




Previous 


Experiment 




LDA 


GGA 


LSDA 


LSDA+SO 




a 


4.6444 


4.7049 


4.6433 


4.6383 


4.6832 6 


b 


2.8505 


2.8946 


2.8467 


2.8437 


2.8717 6 


c 


4.0464 


4.0955 


4.0432 


4.0388 


4.0761 6 


V 


53.57 


55.78 


53.44 


53.27 


54.82 fe 


B 


336.1 


303.45 


364.7 


385.4 


365-395° 



Q see Ref.|33| 
^see Ref.|| 
c see Ref.0 
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TABLE III: Calculated elastic constants of orthorhombic OsB 2 . Note that a relaxation of the 
atomic positions was carried out. All values are in units of GPa. 





Cl2 


Cl3 


C22 


C23 


C33 


C44 


C55 


C&6 


597.0 


198.1 


206.1 


581.2 


142.6 


825.0 


70.1 


212.0 


201.3 
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Figure Captions 

• FIG [1} The ball-and-stick model for structure of orthorhombic OsB 2 . (a) The unit 
cell containing a single formula unit. Two Os atoms at Wyckoff positions of 2a: Osl 
(|, 7, z) and Os2 (|, |, —z). Four B atoms at Wyckoff positions of 4/: Bl (u, j, v), 
B2 (— u + ~, 7, v), B3 ( — u, |, —v), and B4 (u+ |, |, — (b) The supercell repeated 
in three dimensions by 2x2x2, in order to illustate the layers in z(oy c)-direction. Os 
atoms are shown as dark gray balls and boron atoms as light gray balls. 

• FIG 0: Calculated pressure versus fractional unit cell volume for OsB2, along a com- 
parison with available experiment data. 

• FIG |3|: Electronic band structure along the high-symmetry direction of the Brillouin 
zone for OsB 2 . High-symmetry points are labeled as the following ||40|| : T = (0.0, 0.0, 
0.0), Y = (0.0, 0.5, 0.0), S = (0.5, 0.5, 0.0), X = (0.5, 0.0, 0.0), and Z = (0.0, 0.0, 0.5). 
The zero of energy is set as the Fermi level and shown in dotted line. 

• FIG [|: (Color online) Total, site and angular momentum-decomposed density of states 
(DOS) of OsB 2 . 

• FIG [5]: Charge density at the equilibrium lattice parameters of OsB 2 , along the lines 
between (a) Os and B atoms, (b) Os and Os atoms, and (c) B and B atoms. 

• FIG § Contours of the charge density at the equilibrium lattice parameters of OsB 2 , 
in the planes (a) with normal vector along <100> and central point at Osl atom, (b) 
with normal vector along <100> and central point at Bl atom, (c) with normal vector 
along <010> and central point at Osl atom, (d) with normal vector along <001> and 
central point at Osl atom, (e) with normal vector along <001> and central point at 
Bl atom, and (f) through Bl, B2, B3 and B4 atoms. Charge density is in an increment 
of 0.1 e/A 3 from to 2.1 e/A 3 and all distances are in A. 
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